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D . Abstract 

We study the impact on the unitarity triangle of the 5^ - 5^ flavor oscillation data anal- 
^ ■ yses being carried out at the Tevatron Run II. In the current version of this work we con- 

sider the first direct experimental bound on the mass difference Am^^ recently reported 
by D0 along with the first measurement achieved by CDF, Ams^ = (17.33 ±0.4) ps~^. 
In particular we show how these measurements further constrain the compatibility be- 
tween the experimental value of sin 2/9 '^''p- = (0.591,0.783) at 99% C.L. and the value 
obtained by the fit of the CKM parameters considering just CP conserving experimen- 
tal observables: gin 2/3 = (0.742,0.856) at 99% C.L. We also obtain bounds 
on certain classes of processes beyond the SM. 
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1 Introduction 



We are at an exciting era of 6-flavoured hadron physics, as the ongoing experimental 
efforts in this area can be used to test to a better accuracy the CKM matrix and hence 
the Standard Model (SM) giving at the same time bounds on processes beyond the SM 
(BSM). 

Various high precision B physics measurements have been made possible notably 
by the B factories operating at the T(45') resonance. Mixing in the system, with 
the light B mesons being there copiously produced, has been measured accurately. 

The situation is different for the heavier neutral B meson system, B^, which has 
been intensively studied by the Tcvatron experiments, CDF and D0. For the last few 
years they have been accumulating large B^ data samples, developing and optimizing 
complex data analysis techniques which finally lead to the measurement of mixing in 
this system. A task that has been rather challenging for almost two decades. 

All combined experimental data prior to Tevatron Run 11 have given only a lower 
bound on the mass difference that characterizes flavour mixing in the B^ — B^ system, 
namely Atjib^ > 14.5 ps^^. 

The first B^ mixing results from both of Tevatron collaborations after Run II, 
namely Am^^ > 16.6 ps~^ with a combined sensitivity of 20.0 ps~^, were presented 
last year [1] and translated already into sizable contributions to the world average. 
We have studied its effect on the determination of the unitary triangle (UT) in the 
previous version of this work. 

More recently after having extended their AniB, analyses to 1 fb~^ of data, both 
collaborations at Tevatron have reported results of unprecedented relevance. D0 
has established a first direct experimental two-sided bound Am^^ G [17,21] ps^-*^ 
at 90% CL [2]. CDF presented shortly afterward a first and precise measurement 
AmB, = l7.3ltlii ps-i [3]. 

How then these novel results translate not only in the determination of the CKM 
parameters and their compatibility with each experimental measurement, but also in 
constraining models BSM? This is the motivation of the present work. 

The CKM fitter group [6] and the UT fit group [7] have been developing state-of- 
the-art analyses of the CKM parameters which are also regularly updated with new 
experimental information. However we perform this analysis to study the impact of 
AniB, on the determination of the unitary triangle (UT) parameters in the SM and 
also to have a code compatible with such reference analyses in order to test models 
BSM. 

To this end we first update the UT fits of the SM comparing the bound (corre- 
sponding to the combined oscillation data before the CDF and D0 results at Run II 
had become available) and the current measurement. Then we begin our tests of mod- 
els BSM within the Minimal Flavour Violation (MFV) scenario and briefly mention 
consequences for an example of a horizontal symmetry. 

During the last and the present decades many groups have been developing anal- 
yses in order to test the unitary triangle and CP violation using various approaches. 
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including the Bayesian method [8], the Range-fit (Rfit) method [9], based on a frequen- 
tist approach, a Gaussian method [10], and the scan method [11]. These approaches 
treat in somewhat different ways the available information on the experimental and the 
theoretical uncertainties. We use the Bayesian approach to construct a global inference 
function for the Wolfenstein parameters p and f/, from which a posterior probability 
density function (pdf) for all the fitted (within the analysis) and related parameters to 
the unitary triangle and the constraints used to determine it can be derived. 

The analysis of the UT and CP violation tests the SM, hence giving it a more robust 
success or, if physics BSM is present already at sizable levels, detecting inconsistencies 
with the various unitary tests. In what it has been called the Classic analysis of the 
UT [40], the CP conserving parameters - |Ktfe|/|Vc6|, ^ttlb^ and AmB^/AmB^ - and 
CP violating ones - ck and sin 2/3 - are taken into account. The ratio iKibl/lKbl 
is not expected to have a big contribution from beyond tree-level processes of semi- 
leptonic decays from which ]T4&| and \Vcb\ are extracted. The CP asymmetry in S — > 
J/ipKs could give deviations from measuring not just the SM unitary angle [3 but 
P + Obsm- Over the past years the SM prediction of this quantity using the other 
UT parameters was consistent with its observed value. But an important result of the 
present analysis (as well as other recent analyses) is that there is an incompatibility 
with the input experimental value sin 2/3 — (0.623, 0.751) at 95% confidence level 
(CL) and the output value of the fit including just the CP conserving parameters: 
sin 2(3 cp-conscrv. ^ (0.760, 0.836) at 95% CL. Now they are just compatible at 99% CL.^ 
This is an important hint for BSM processes but of course a better measurement of 
Atub^ and a better determination of the parameters \Vub\ and \Vcb\ will be needed in 
order to shed light on this apparent incompatibility. 

We just consider AF = 2 fiavour changing processes, leaving out of our analysis 
AF = 1 processes, which despite being relevant in some cases the present experimental 
bounds are not at the level of precision [14] as most quantities entering in the Classic 
analysis. 

Mixing in the kaon and in the Bd,s neutral meson systems, from which €k, ^f^Ba 
and Ams^ are respectively obtained, can be particularly sensitive to physics BSM. 
The most general way to identify the deviation of e/^, Am^^ and Am_B^ from their SM 
prediction is by measuring the deviation from unity of the ratios 

Im{(X°|7Yj^*^|X°)} / Im{(X°|:^^SM|^o^| 

respectively, where Ti.Jg''^^ is the effective Hamiltonian that takes into account all the 
processes, both in the SM and BSM, of the relevant AF = 2 (X°-X° and 5° ^-5° J 
flavour changing neutral current processes. 

However, we have to bear in mind that there are two main restrictions on measuring 
these ratios. The first one is that there are still big uncertainties in some quantities 

^Using the bound Eq. (20) they were compatible at 95% CL output since the value of this case for 
the CP conserving processes is sin 2/3 = 0.793^o.024- 
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entering in the computation of the effective SM Hamihonian matrix elements, chiefly 
Bk and Jq^Bb^, which obscure differentiation of these and the effects of processes 
BSM. The second one is that in general these ratios are related in different extensions 
of models BSM and just for the minimal flavour violation (MFV) scenario of the Min- 
imal Supcrsymmctric Standard Model (MSSM) the form of the effective Hamiltonians 
and the QCD next to leading order (NLO) determination of its parameters has been 
computed [31]. 

For the flrst of the aforementioned restrictions not much can be done currently but 
for the second one we can test models for which we can extract meaningful information 
from the current way in which the UT analysis it is performed. Since we want to test 
what is the impact of the new results on Auib^, we explore the cases of physics BSM 
for which the ratio Atub^/ AruB, is equal to the case of the SM. 

The MFV violation scenario, which essentially requires that all flavour and CP vio- 
lating interactions are linked to the known Yukawa couplings, is a well known example 
of this condition. The MFV scenario can be easily implemented in the MSSM, as it has 
been done extensively since the pioneering work of [15], but using an effective operator 
approach [16] it can also be implemented for other cases. With the current bound 
we can see that the MFV supersymmetric scenario appears plausible with the current 
Attibs value. 

Given the new experimental information we thus study the possibility of having 
significant contributions of processes BSM in the determination of the CKM elements, 
and we show that, under certain constraints, somewhat large effects of processes beyond 
the SM are still possible. We mention briefly the consequences of this analysis for 
models of fermion masses with an SU{3) horizontal (flavour) symmetry. 

2 Constraints of the CKM matrix 

Here we just quickly state the formulas used in the SM fit of the CKM matrix, V, 
and remark the processes from which they are extracted in order to point out which 
constraints we use to put bounds on processes beyond the SM. The inputs of the fits 
are presented in Tables (2)-(3) of Appendix A. We recall here that V can be expressed 
in terms of the four Wolfenstein parameters: A, X p and rj. All elements of V can be 
obtained as an expansion in A up to the desired n-th order, using the unitary condition 
from the basic definitions 

IK.I = A, \Vcb\ = AA^ iKbl = AX'ip - irj), (1) 

which are valid to all orders in A. The parameters A and A are obtained from tree-level 
kaon and bottom semi-leptonic decays which have large branching ratios making their 
determination independent, to an excellent approximation, of any processes BSM. The 
re-scaled parameters p = cp and fj = crj, with c = 1 — A^/2, arc the coordinates of the 
apex of the unitary triangle (UT), with sides Rh and Rt defined by 

Rl = p' + f, Rl = (1 - pf + f (2) 
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and angles a, (3 and 7 defined by 



sin/3 =—, sin7=— , a = 7i-p-j. (3) 

Kt Kb 

In the Classic fit of the SM CKM matrix, p and fj are extracted from 
(i) CP conserving processes: & (ii) CP conserving processes: 

j^, Ams, and and sin 2/3. 

In (i), |T4;,| and \Vcb\ are extracted from semi-leptonic decays of B mesons. The pa- 
rameters ArriB^ and ArriB^ arc extracted from measurements of B^^-B^^ oscillations, 
which arc = 2 flavour changing neutral current processes (FCNC). From (ii), the 
asymmetry parameter |e/^| is induced by K^-K^ transitions, which are A^* = 2 FCNC 
processes. Finally the measurement of sin 2/3 is dominated by the channel B^ — > J/i/jKs 
and complemented by other decay modes. 

The parameters {exl, \Vuh/Vcb\, Anis^, ArriB^ and sin 2/3 are expressed directly or 
indirectly in terms of p and f], and a set of theoretical and experimental parameters x, 
some of which have been well measured and some have still large uncertainties. Thus, 
in the fits the parameters x which have been well measured are kept fixed and those 
with dominant uncertainties are varied and fitted. 



\Vub\/\Vcb\ 

The CKM matrix elements \Vub\ and \Vcb\ are measured in exclusive and inclusive 
semi-leptonic B decays. In terms of the re-scaled Wolfenstein parameters, |T4&|/|Kb| is 
expressed as: 



\Vcb\ c 



(4) 



Within the SM, the mass differences Arud^g of the oscillating systems B^^-B^^ are 
very well approximated by the relevant electroweak transition ("box") diagrams, de- 
scribed by the Inami-Lim functions [17], which are dominated by top-quark exchange. 
Am^rf^ arc proportional to the moduli of the matrix elements of the effective Hamil- 
tonian for the AB — 2 transitions, 

^^B, = ^,rnBjlBs,VBSix,)\V*Vj' oc \{B,\n!^\B,)\, {q = d,s) (5) 

where Gp is the Fermi constant, rjB is a QCD correction factor calculated in NLO, tub 
and rriw are the Bg meson and W boson masses respectively. 
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The dominant uncertainties in Eq. (5) come from the evaluation of the hadronic 
quantities: the B meson decay constants Jb^ and the bag parameters Bb^, which 
parameterize the value of the hadronic matrix element. The Inami-Lim function, S{xt) , 
describes the | AS| = 2 transition amplitude in the absence of strong interaction, where 

2 

the mass of the top quark enters via Xt = jt^- ^ 

In terms of the Wolfenstein parameters, Eq. (5) above for the B^ case can be 
expressed as 



Ams, = C^^A'\%1 - -pf + rf]mBjlBB,VBS{xt), 



(6) 



where wc have defined the constant Ca™ = q^2^ ■ The parameters with dominant un- 
certainties in Eq. (6) are Jq^Bb^, A and A; among these, A and A are varied parameters 
of the fit. 

The size of side | Vtd|/(A| Vcbl) of the unitary triangle can be obtained from the ratio 
of Amfi^ and Am^^, 



Am 



Am; 



ruB, 


Vtd 


2 


itibA^ 


Vts 


2 



with ^ 



(7) 



which is expected to be less dependent on the absolute values of the non-perturbative 
quantities /s and Sb, and hence it can be characterized by ^, an SU(3) breaking 
correction whose value is obtained from lattice QCD calculations with better precision 
than the individual matrix element quantities themselves. In the fit we use ^ and 



fssV^Bs instead of fB^\J Bb^^, which makes the constraint Am^^ more effective [41]. 
The constraint we use from AniB^ is expressed, from Eq. (7), as 



AmR = Amfi, 



rriBs 



rriBa A2 (1 - p)2 + r] 



^2 ' 



(8) 



where Am^^ is here taken as an experimental input. The parameters with dominant 
uncertainties in Eq. (8) are ^, A, A, which are varied parameters of the fit. 



The parameter ex expresses the measurement of indirect CP violation in the neutral 
kaon system. In terms of the Wolfenstein parameters it is given by 



Im Asm{K) 



C^BkI^ Asu{K) cx {K^\Hl^'\K^) , 
A''X%-7],cXc + A^X^ (1 - p) r]uS{xt 



with 

+ r]ctS{xc,Xt)] , 



(9) 
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The short distance QCD corrections are codified in the coefficients rjcc- Vtt and rjct, 
and are functions of the charm and top quark masses and of the QCD scale parameter 
Aqcd] these parameters have been calculated in NLO. The Inami-Lim functions, which 
describe the \A.S\ — 2 transition amplitude in the absence of strong interactions, are 
given by S{xt), as in the Am^^ case ^, and by S{X(,. Xt). ^ 

The parameters with dominant uncertainties are Bk, rjcc, Vet, 'f^c and rrit which are 
hence varied in the fit. 

sin 2^5 



A direct determination of the angles of the unitary triangle can be achieved via 
measurements of CP asymmetries in various B decays. These angles are defined in 
reference to the normalized unitary triangle: 



(10) 



a — arg 


[ VudV:, \ 


, /3 = arg 


\-VcdV:A 


, 7 = arg 


[ VcdV;, \ 


The value of 


sin 2/3 is 


measured in 


h CCS transitions, 


such as 



decays. The UT angles can be expressed directly in terms of the re-scaled Wolfenstein 
parameters; in particular we have 



sin 2/3 = 



2^(1 -p) 
+ (1 - p)2 



(11) 



How to account for indications of physics beyond 
the SM 



The ratio iKbl/lKbl is not expected to have a big contribution from beyond tree- level 
processes of scmi-leptonic decays from which \Vub\ and \Vcb\ are extracted. Mixing in 
the kaon and B^^s systems however can be particularly sensitive to physics beyond the 
SM (BSM). Their deviation from the SM can be expressed, without loss of generality, 

by 





-T/Total 


\Bd,s) 






Bd,s) 



lm{{K' 



1 -^Total I 1^0 



'-cff 



lm{{KO\n^^\KO)} 



(12) 



New CP phases may arise from the interference between mixing and decay ampli- 
tudes in the Bd,s systems. Thus CP violating asymmetries would involve not only a, 
(5 and 7 but additionally Ob^i and Ob^i defined by 



sin 2/3 = sin(2/3^^ + 2^B. 



a 



a 



SM a 
— U} 



(13) 



S{Xc, Xt) — Xc In Xc -h Xc 



4(l-xt)^ ^^^^ -f- 4 



xt—i 
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and similarly for the B'^ system whose phase information may be obtained from mea- 
surements of CP asymmetries in — > J/ilxj) which however are not yet available. For 
the kaon system, ArriK is not considered due to the lack of control in the QCD long 
distance effects in the K'^-K^ system and hence a possible additional phase BSM in 
this sector is neglected. In general the BSM contributions to the amplitudes of the 
processes responsible for Auib^ and sk are no longer proportional to the CKM matrix 
elements as in the SM, 

{B,\n!^\B,) oc {Vt,V*r, iq = d,s), 

{K'lntm') oc {V.aV:,r,{VuVZr,{V^V:,){VuVZ), (14) 

but at best there will be only additional contributions with the same phase. 

This is the point that can be used if we want to extract useful information of 
processes BSM from the analysis of the unitary triangle. Meaningful bounds on such 
processes can be extracted from B^^-B^^ and K^-K^ as long as the corresponding 
effective Hamiltonians are still proportional to the combination of CKM elements as 
given in Eq. (14), and such that there is only a change in the contribution to the box 
diagrams of these AF = 2 transitions. This amounts to a replacement of r)2S{xt) in ex, 
as appears in Eq. (9), and a replacement of riBS{xt) in AmB^ and ArriB^, as appears 
in Eq. (5), by new functions describing the contribution of particles BSM. These can 
be generically written as 

ArriB, = -^^mBjlBBM*M'rjB[S{xt) + SS^rnsJ, and 

le^l = CBkA^X^t] [-ri^x, + rj^tSix,, Xt) + A^\\l - p)riu[S{xt) + 6S,j] (15) 

respectively, by normalizing the QCD correction factors of the processes BSM with 
respect to those of the SM. 

The set of theories which satisfy these conditions is quite restrictive, but it is im- 
portant to note that the bounds put on them are at the same level of precision at which 
the SM unitary triangle parameters can be currently determined. On the other hand 
we cannot test by these means new exciting theories for which there are contributions 
to box and penguin diagrams that are not proportional to the CKM matrix elements as 
the SM top contribution, or there are complex phases beyond the phases of the CKM 
or new local operators (contributing to the relevant amplitudes) and hence introducing 
additional non-perturbative factors B^ and box and penguin diagrams. 



4 Statistical analysis of the unitary triangle 
4.1 Likelihood technique 

We briefly introduce here the general hues of the inference method we adopt, mostly 
for completeness and for reference in later sections. The general approach is well 
documented in [18] where further generic clarifications may be found. 
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The likelihood technique for the analysis of the unitary triangle in the SM has 
been widely studied. The Bayesian approach has been investigated by the UT Fitter 
group [7] . We have chosen to follow such an approach in implementing our code, while 
placing emphasis on the importance of the increasing knowledge of Am^^. 

In the Bayesian approach, the combined probability distribution for p and 77 is 
identified with the likelihood 



where f{cj\cj{p,fj,{xi})) is the conditional probability density function (pdf) of the 
constraints Cj e {|146|/|K6|) Ams^, AniB^/ ^tub^, e^r, sin 2/3}, given their dependence 
as functions of the CKM parameters, p, f], A, A and Xi. Besides the constraints 
themselves, there are two classes of parameters involved: (i) fitted, for which pdfs 
fi{xi) are constructed (e.g. the top mass); and (ii) fixed, which are taken as constant 
(e.g. the W mass). 

The jjnor probability density functions f{cj) for the constraints are considered to be 
Gaussian distributions, possibly convoluted with flat distributions. The exception was 
AniBg, prior to its first measurement, for which a more complete set of experimental 
information was used. 

More specifically, the pdf used in this case, in order to be consistent with the 
Bayesian approach, is implemented via the likelihood ratio, R Eq. (23), after accessing 
the amplitude point (A, a^) associated to the frequency value obtained by evaluating 
the r.h.s. of Eq. (8). The method is described in further detail in a following section. 

For fi{xi) we consider Gaussian and in some cases flat probabilities. A Gaussian 
pdf is chosen when the uncertainty is dominated by statistical effects, or there are 
many contributions to the systematics error, so that the central limit theorem applies. 
Otherwise, a flat (uniform) distribution is used for the uncertainty. When both Gaus- 
sian and flat uncertainty components are available for a parameter, the resulting pdf 
is obtained by convoluting the two distributions. I.e., for an observable parameter x of 
true value x, with Gaussian and uniform uncertainty components, ag, au, one has for 
the parameter and its pdf, f{x), 

X = X + Xg + Xu f{x) = 5{x — x)^ Gaus(x|(jg) (8) Unif (x|(jf). 

The integration of Eq. (16) can be performed using Monte Carlo methods, followed 
by its normalization. Hence we can calculate the resulting pdfs for parameters involved 
in the flt. 

The expression of Eq. (16) shows exphcitly that whereas a priori all values of p and 
77 are equally likely by assumption, i.e. fo{p,v) — const., a posteriori the probability 
clusters in a region of maximal likelihood. 

The probability regions in the {p,f]) plane are constructed from the pdf obtained in 
Eq. (16). These are called highest posterior density regions, and are defined such that 
C{p,fi) is higher everywhere inside the region than outside. 




n /(ci|ci(p,r7,{a;i})) X JJ fi{xi) dxi x fo{p,r)). 



■1,M i=l,N 



(16) 
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:= {z = {p,r]) : Jp^ C{z)dz = w; C{z') < minp^C{z),W^,^p^}. 

The one dimensional pdfs are obtained similarly. For example, the p pdf is obtained 
as C{p) oc / C{p,fi)df], from which its expected value can be calculated together with 
the corresponding highest posterior density intervals. 

A similar procedure can be in principle used in order to obtain the pdf for other 
desired parameters. Technically, one may also use the probability function for trans- 
formed variables; i.e., that for u(x) one has /(u) = /(x)|9x/9u|, where the last factor 
denotes the Jacobian. In this way, the pdf for a parameter x can be obtained from 



where C{p,f]) has been computed in Eq. (16) above. 

Besides the probability distribution in the (p, 77) plane, we are most interested in 
obtaining the posterior probability distribution for the Athb^ parameter. 

5 Using the available Bg flavour oscillation results 
5.1 Experiments — short overview 

The study of flavour oscillations in the B^-B^ system has been the subject of many 
past experimental analyses, performed at ALEPH, CDF, DELPHI, OPAL and SLD. A 
wide range of data analysis techniques has been employed, which have been developed 
for taking advantage of detector capabilities and the characteristics of the collected 
data samples. These include from fully exclusive, where the Bg decays are completely 
reconstructed, to fully inclusive methods, which aim at identifying the mesons decay 
vertices; along with several 6-flavour tagging techniques. The most sensitive among 
these analyses were those based on semi-exclusive or semi-inclusive lepton samples, 
where the Dg, resulting from the semi-leptonic Bg decay, was either inclusively or 
exclusively reconstructed. 

In general, the more inclusive analyses benefit from considerably larger number of 
candidates, while the more exclusive ones take advantage of the more precise infor- 
mation about the B decay candidates which is available. In practice, therefore, the 
resulting sensitivity to B oscillations is dependent on a trade-off between the quantity 
and the quality of the events forming the samples. Furthermore, the relative weight 
of the latter increases quickly when probing higher oscillation frequencies. Accord- 
ingly, the significance of an oscillation signal at relatively lower frequencies benefits 
more readily from high statistics samples, provided the events are characterized by 
sufficiently adequate resolutions. This has to some extent been the case until rather 
recently. The situation becomes different for a signal of a higher frequency, as is known 
to be the case of the B^ system — a measurement of Atub^ requires samples with very 
good effective resolutions together with adequately large yields. 

A new generation of B^ mixing analyses is being undertaken by the CDF and D0 
collaborations, which are accumulating large samples of B^ meson decays produced 




(17) 
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in the pp collisions of the Tcvatron at Fermilab. Both collaborations reported first 
preliminary results during the year of 2005 [1]. 

Samples of partially reconstructed semi-leptonic Bg — ^ Dglu decays have been col- 
lected by both experiments, whereas CDF has also gathered substantial samples of 
fully reconstructed Bg — > Ds{'K'k)'k decays. The Dg meson is reconstructed in several 
exclusive channels. 

The partially reconstructed decay samples have originally provided the predominant 
contributions to the combined sensitivity, in view of the considerably larger yields. 
Nevertheless, as the sizes of the fully reconstructed decay samples increased, these 
have started to provide unprecedented effective resolution and large sensitivities at 
higher probe frequencies. 

The preliminary results reported last year have contributed already significantly 
to the world average [12], increasing the combined sensitivity to 20 ps~^ and further 
pushing the lower exclusion limit. 

Tevatron projections based on the those preliminary analyses have demonstrated 
the potential in Run II for covering the SM favored region for AuiBs- Indeed this 
has already taken place most recently as, after using optimized analysis techniques 
and increased datasets, sensitivities close to 30 ps~^ have been reached. And indeed 
evidence of oscillation signal within the SM favored region has been estabhshed. 

Following the Tevatron the excitement of the study of the B^ system will be trans- 
ferred to the LHC experiments: the LHCb, as well as the general purpose detectors 
ATLAS and CMS. A confirmation of the Am^^ measurement, with an increased preci- 
sion, is readily expected. Additional analyses will be pursued, which will provide crucial 
information for instance on the phase of the mixing amplitude and on CP violating 
processes. 

5.2 Measurements — amplitude scan 

The study of oscillations is based on the analysis of the proper decay time t of flavor 
tagged Bg candidates. 

It involves the measurement of the B meson decay distance and momentum, as well 
as a determination of whether the particle-antiparticle system decayed with the same 
or the opposite flavor with which it was originally produced. 

The model used to describe the proper time distribution involves the exponential 
decay dictated by the system's lifetime r, modulated by an oscillating term describing 
the probability for the Bg to have mixed, 

P{t) ~ - exp^r (1 - AcosAuibJ) ■ (18) 

T 

This expression needs to be adapted to take into account detector resolutions, recon- 
struction effects, and imperfections of the &-flavour tagging methods used. 

The parameter A has been introduced in 18 in an ad-hoc fashion for accessing the 
amplitude of the oscillations using the so-called amplitude method [13]. Instead of 
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extracting from the data the quantity of interest - Am^^ - directly, in this method 
the hkehhood is maximized as a function of A, while Am^^ is a fixed parameter of 
the fit. The procedure is repeated for many different probe values of Am^^, and the 
result is thus conveyed as a set of measured ^(AmeJ and (T^(AmBj in a spectrum 
of probed oscillation frequencies. The value of A should fluctuate around zero unless 
the true oscillation frequency of the system is being probed, for which case A should 
be consistent with unity, within errors. 

A given frequency value is excluded at 95% CL, following a one-sided Gaussian test, 
if the corresponding amplitude and its uncertainty satisfy A + 1.645(7_4 < 1; the lower 
exclusion limit of an experiment is defined as the value below which all frequencies 
satisfy the latter condition. The sensitivity of an analysis is defined as the largest 
frequency for which 1.645(7^ < 1 holds; i.e. it is the largest excluded frequency if A 
were exactly zero as it is expected to be the case in the absence of a signal, removing 
the dependence on fluctuations of the central value A which affects limits. 

The amplitude method was originally proposed [13] as a convenient approach to set 
limits, as well as for combining such exclusion regions obtained with different analyses. 
It has provided a consistent fitting procedure through which all experiments have com- 
monly expressed their results. The combination of amplitude results is standardly done 
most straightforwardly as averages of Gaussian measurements. This is expected to hold 
in general to good approximation; more detailed combinations could be performed for 
instance by combining the likelihood profiles for each A measurement, which however 
would not be feasible as it would be based on detailed information which is no longer 
available for previous measurements. 

In the averages performed by HFAG [12] the measurements are adjusted on the 
basis of common physics input values, and possible statistical correlations are taken 
into account. 

Prior to the inclusion of the Tevatron Run II results, the world-average exclusion 
limit was [12] 

AruB, > 14.5 ps"^ (95% CL) , (19) 

with a sensitivity of 18.2 ps~^. 

With the inclusion of the Tevatron results reported in the past year (2005) [l],the 
limit was increased to 

ArriB, > 16.6 ps"^ (95% CL) , (20) 

with a combined sensitivity of 20.0 ps~^. 

More recently, after the analysis of the first fb~^ of Tevatron data, both collab- 
orations have reported their corresponding results. D0 has published a first direct 
double-sided experimental bound [2] placing Am^^ in the interval 

17 < AruB, < 21 ps"^ (90% CL) ; (21) 
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the measured sample sensitivity is 14.1 ps ^. CDF has pubhshed a first direct, precise 
measurement [3], 

ArriB, = n.Sltofsistat.) ± 0.07{syst.) ps"^ , (22) 

having been extracted from data samples whose measured sensitivity is 25.8 ps^^; 
the measurement has aheady a 2% level precision. The probability for the observed 
signatures being due to background fluctuations were evaluated by the collaborations 
to be of the order of 5% (D0 ) and 0.2% (CDF). This is still below the threshold 
commonly set for attaining the status of observation (Sa); a (potential) confirmation 
of the measurement with an even larger significance is being awaited. 



5.3 Likelihood implementation 

In the fit to the CKM parameters we use the complete information about Am^^ which 
is provided by the full amplitude scan. 

The measured values of the amplitude and its uncertainty, A and a^, may be used to 
derive [13], in the Gaussian approximation, the log-likelihood function, Aln£°°(AmBj, 
referenced to its value for an infinite oscillation frequency 

AlnC^iAruBj = In £(00) - ln/:(AmBj = (I - A) \, 

V2 J CT_4 

Alni:°"(AmBj^i, = Aln/:°°(AmsJnomix = ■ 

The last two relations give the expected average log-likelihood value for the cases when 
ArriB^ corresponds to the true {mixing case) or is far from {no-mixing case) oscilla- 
tion frequency of the system, characterized respectively by unity and null expected 
amplitude values. These are shown in Figure (1 i). 
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The log-likelihood difference, according to the central limit theorem of likelihood 
theory is x^-distributed: AlnC = |x^. We translate therefore the amplitude scan into 
the likelihood ratio 



£(oo) 




(23) 



through which the constraint for Ame, is implemented in the fit. This is represented 
in Figure (1-ii). 

Wc rc-itcratc that the exponent in Eq. (23) corresponds to the x^, or log-likelihood, 
difference between the cases where an oscillation signal is present and absent, for which 
the true amplitude value is 1 and 0, respectively. Note that hypotheses for Athb^ 
associated with larger and more precisely measured ^-values in the scan contribute a 
larger weight in the fit. 

Extending the amplitude spectrum 

The amplitude fits are usually performed for Auib^ values lower than a given thresh- 
old, beyond which the fit behavior may become unstable. On the other hand, in the 
framework of the CKM fit, it is in principle desirable to have R defined for all positive 
frequency values, which in turn demands for a continuation of the amplitude spectrum 
beyond what is measured. 

The extrapolation of the value of may be achieved through an analytical, Fourier- 
based description [13] of the measured significance curve. 



where at denotes the uncertainty in the measured proper decay time (this expression 
would in fact need to be properly averaged over the samples signal uncertainty distri- 
butions) . 

An extrapolation of A itself as such is not possible, however it is here sufficient 
to note that its expected values lie in the vicinity of either zero or unity as already 
mentioned. Therefore, for the unmeasured part of the spectrum, the exponent in 
Eq. (23) becomes small, and quickly approaches zero. We take this asymptotic limit 
as the criterion for extending R beyond the experimentally probed frequencies. 
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Figure 2: Probability density function of AniB^ as a result of the unitary triangle fit using 
(i) the data corresponding to Aijib^ > 16.6 ps bound, and (ii) the data corresponding to 
{AtuBs = 17.33 lb 0.4) ps (iii) shows the output of Am^^ without using its constraint in 
the fit. 
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6 Tests of the SM in the UT analysis 

6.1 Classic fit with bound and /irst measurement of AmB 




(ii) 

Figure 3: The first (i) of these plots represents the (/), ij) plane with the bound results and 
the second one (ii) with the data of the first measurement on Am ■ 

Although AF = 1 processes such as the penguin transitions occurring in the charm- 
less decays B —* tttt, pn, pp, which determine a, and the CP asymmetries in semi- 
leptonic B decays {B l~X) are important, we do not consider these constraints here. 
The experimental constraint on a through the mentioned decays is currently sufficient 
even without necessarily assuming the SM, but the relation a = 7T — (3 — jis still used 
for this bound where 7 here comes from the decay amplitudes and (3 from the B'^ — 5° 
mixing, so this AF = 2 process imposes already a good constraint on a. 

The CP asymmetry in semi-leptonic decays, Asl, is a crucial constraint of the UT 
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analysis because once r^^ and 9b^ arc defined through Eq. (12), the generahzation of 
AsL to account for BSM processes depends on r^^ and 9b^ [40] at the NLO in the 
penguin term in the AF = 1 amphtude. However the present experimental bounds are 
not precise enough to put constraints on p and f] [14]. 

The two main fits that we have performed include the constraints |Vu6|/|V^6|, Am^^, 
AniB^/ AuiB,, €k and sin 2/?. The first one (i) uses the bound Eq. (20), and the second 
one (ii) uses the first measurement Eq. (22) by CDF.^ 

In Figure (3)(i) we present the 68%, 95% and 99% CL of the two dimensional 
pdf for p and fj with the values of Tables (2) and (3) using the bound Eq. (20) and 
in Figure (3)(ii) using the current measurement Eq. (22). 

As we can see in the (p, f/) plane, the radius of the circle describing the constraint 
Ams^/Ams^ has been reduced and with it the overlap region of all the constraints, 
specially for iVubl/lKftl and sin 2/3. As a consequence the central value of p has been 
increased and the central value of fj has been shghtly reduced. In Figure (4) and 
Figure (5) we present the pdfs of p and fj for these cases. As we can see by comparing 
these figures, the one dimensional (Id) pdf of p has been shifted to the right while the 
pdf of fj has been left practically unchanged. In Figure (2) we present the Atub^ Id 
projection using the (i) bound, (ii) the measurement, and (iii) the fit output of Am^^ 
without using this constraint. 

As expected with the current measurement of Am^^, the overlap region within 
the allowed 68% CL can determine with a smaller uncertainty the values of p and fj. 
The last bound of Auib^ was closer to the region of maximal likelihood of the other 
constraints, and the relative uncertainties of p and fi were 15% and 9% respectively, in 
contrast to the current 11% and 4%. 

We also compare how well the parameters defining a unitary triangle, x — Rb, Rt, a, 
P and 7, are fitted in each of the cases (i) and (ii), by a relative to the experimental 
value, 

2 



where x is computed as given in Eq. (2) with the values of p and fj returned by the fit, 
as shown in Figures (4, 5), and x represents the values as shown in Figures (6) and (7). 
We find that the fit including the measurement of Am^^ (ii) is fitted slightly worse to 
the constraints on unitary than the bound (i) fit: 

Xl,, = 2.3xl0-^ xl,, = 5.9xl0-\ 
Xi« = 0.008 Xi,, = 0.018 

Xsin2/3(,) = 1-1 X 10-^ Xsin2/3(«) = 1-4 X 10"^ 

^In the original version of this work we used the available information from the s'^d s oscillation 
experiments for the summer conferences of last year (Am^^ > 14.5 ps~^) and other fit using the data 
after the summer and fall conferences (Atob^ > 16.6 ps^^), averaged by Heavy Flavor Averaging 
Group (HFAG), using the information by the CDF, D0 , ALEPH, DELPHI and OPAL coUaborations. 
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X; 



sm7{i) 



1.2 X 10" 



5.8 X 10"^ 



(25) 



Since a is determined through the constraint a = vr — /? — 7, the difference in sin a has 
no meaning for this kind of comparison. 




0.1 0:2 0:3 0:4 




0:2 0:3 0.4 0:5 0:6 



Figure 4: p and r/ one dimensional projections using the bound results Eq. (20). 





posterior for T[ 
68 % C.L. interval 
T = 0.363;r. 



P 1 

Figure 5: p and f/ one dimensional projections using the results of the fit with the Am^^ 
measurement Eq. (22). 
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Figure 7: Rh and Rt one dimensional projections using the data of the Atub^ measurement. 




Figure 8: a and 7 one dimensional projections using the bound results. 
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sin(2|3) sin(2p) 



(i) (ii) 
Figure 10: sin 2f3 one dimensional projections using the (i) bound and (ii) the fit using 

the data of the fit with the Amn measurement. 



20 



6.2 Indirect evidence of CP violation 



The analysis of the unitary triangle (UT) allows for the comparison of the fit using 
only the CP conserving processes, |Kife|/|Kb|, ^t^B^ and Ams^, with the current ex- 
perimental values of the CP violation parameters ex and sin 2(3. This is a way to check 
how consistent is the fit of the parameters that are sensitive just to the sides of the 
UT with respect to the measurements that are sensitive to the angles, which it is an 
indirect measurement of the amount of CP violation. The (p, f/) plane of this fit is 
shown in Figure (6.2). We also compare the fitted Id pdfs of the parameters {exl and 




Figure 11: {p,fj) plane with the results of the fit using the AthBs nieasurement and using 
only the CP conserving constraints |Ki6|/|V"cfe|, Am^^ and AniB^/ Ams^- In dotted lines we 
have plotted the CP violating constraints: sin 2/? and ex to compare them with the pdf of 
this fit. 

sin 2/5 to those given for the case when these parameters have been included in the 
fit as constraints. The output value of e^- for the fit including just Attib^ 
and AniB^/ AniB^ is ex = (2.230,2.96) x 10~^ at 68% CL; it is consistent with its 
experimental value of ex = (2.11,2.45) x 10"^ at 68% CL. However as we can see by 
comparing Figure (12 a) to Figure (10), for both the Am^^ bound and measurement 
cases, the value of sin 2/3 increases if the CP violating constraints sin 2/? and ex are not 
taken into account: 

sin 2/? CP-conserv. _ ^^g^ q g^g) ^ (26) 

which does leave out the 68% CL interval of the experimentally measured value: 

sin 2/3 '^"P- = (0.655,0.719). (27) 

This inconsistency is already apparent in the fits which include sin 2/? and ex as a 
constraint, but here we note that the values of Eqs.(26) and (27) are consistent just at 
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(i) (ii) 
Figure 12: One dimensional projections of (i) sin 2/3 and (ii) ex using data of the fit including 

the AniBg measurement with out imposing them as a constraint. 



99% CL. At 99% CL, the CP-conserving interval is (0.742, 0.856) while the experimental 
one is (0.591,0.783). 

One reason for observing this apparent inconsistency may be due to the current 
set of experimental inputs for the constraints and their errors, especially for iKib/Kbl- 
A better determination of these semi-leptonic parameters is crucial for checking this 
inconsistency. Another reason is a potential non-zero phase accounting for processes 
BSM in sin 2/? such that we have sin 2/3 = sm{2P^^'^ + 29bj^), as we have discussed in 
Section (3) and as we will be presenting in Section (7.3). 
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6.3 Determination of fg^BB^ and Bk 

The actual measurement of Aitlb^ also provides, within the context of the SM, a 
strong constraint on the non-perturbative QCD parameter f^^Bs,, and through the 
relation Eq. (8) on the parameter Bb^, which suffers from large uncertainties. In 
particular here we can check how its lattice computation, 

II^BbT = (0-223 ± 0.033 ± 0.012) GeV^ (28) 
agrees with the current Atub^ results. 




Figure 13: The pdfs for /| Bb^ [GeV^] and Bk for the bound fit. 




Figure 14: The pdfs for f'^^&B^ [GeV^] and Bk for the fit using the measured Atub^- 

We can read the values of ]3b^ from Figures (13) and (14). The values of this 
parameter (in GeV^) for the overall fit are 

f2 Abound n r)r)Q+0.026 f2 Amoas. n 01 O+0.007 

JBa-^Bd ~ '-'•^^^-0.028' JBa-^Ba ~ '-'•^-'-°-0.008 ' 

= 2.9x10-, xJi^^,^,, = 8.1x10-, (29) 

hence checking that the QCD Lattice value, Eq. (28), and the fit value of Bb^ agree 
remarkably well. 
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Comparing Figures (13) and (14) we can see that the measured 5° — oscilla- 
tion data has had an important impact on /g Bb^, by decreasing its value and its 
uncertainty. This uncertainty is determined to a better accuracy than the QCD de- 
termination, Eq. (28). The impact of Aitlb^ on Bk is not as strong as the impact on 
Jb^Bb^ since it is just affected through the overall fit and not through a constraint, Bk 
only appears in the SM due to the je^^l constraint, Eq. (9). Although the difference, 
SBk, between Bk as obtained in the fit. Figures (13) and (14), and Bk as given by 
the QCD lattice calculations is smaller in case (ii) of measured Atjib/- 

^QCD ^ 0.79 ±0.04, 

Xk,, = 2.26X10-S 

xl,,, = 7.66 xlO-^ (30) 

the uncertainty has increased for the measured Am^^ case. It is interesting to compare 
the effect on Bk when |e/^| is removed as a constraint, as we can see in Figure (15), 
the value of Bk in this case is 




Figure 15: The pdf of Bk when gk is not imposed as a constraint. 

Bk = 0.707t°;°el • (31) 

When the Am^^ bound Eq. (20) was presented, we obtained Bk = 0.745to;o72 which 
agreed quite well with its QCD value, Eq. (30), although the error was of 12% in 
contrast with the 10% error of the QCD caculation. Now we see that although the 
error has decreased to 10% and there is compatibility with the QCD calculation at 
68% CL, here we have Xbj^ = 4-3. 

7 Bounds on processes BSM 
7.1 Scenario I 

In this scenario we consider the constraints on models for which the conditions of what 
has been called Minimal Flavour Violation (MFV) are valid. An effective theory is said 
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to be minimal flavour violating if all higher- dimensional operators constructed from the 
SM and fields responsible for giving the corresponding structure of Yukawa couplings 
(e.g. flavon -scalar- fields whose vevs determine the effective Yukawa couplings) are 
invariant under the CP symmetry and the fiavour symmetry, determined by the group 
Gt. When this happens the dynamics of fiavour violation is completely determined by 
the structure of the SM Yukawa couplings. If all Yukawa couplings are small except for 
the top, then the only relevant non-diagonal operators entering in the Hamiltonians 
of the AF — 2 transitions are proportional to [YyY^]. ~ Ut^u^tj- By looking at 
Eq. (14) we can see that VtdV^ and V^^^V^*, the main contribution to the mixing in 
the AS = 2 and /S.S — 2 transitions respectively, are multiplied by the same Inami- 
Lim function S{xt) (as expressed in Eq. (5) and Eq. (9)) exactly because of the top 
dominance). Once there is another particle contributing at the same order of the 
top contributions then the expressions of the dominating diagrams for the processes 
AS = 2 and AS — 2 become different. This is evident in supersymmetric theories 
when tan /3s is large because the bottom contribution becomes relevant and hence it 
gives a sizable contribution to the AB — 2 transitions. Due to the potential difference 
in these processes, two possibilities are considered: 



7.1.1 Example: supersymmetric MFV scenario in the low tan /3s limit 

In the context of MSSM models there are four physical phases: ^ckm, = arg, 6^ ~ 
arg(//) (where A is the common trilinear term between left and right s-fermions and a 
Higgs field and ji the bihnear term between the two MSSM Higgs fields) and the QCD 
vacuum parameter 9qcd which however can be conveniently set to be zero. 

Experimental upper limits on the electric dipole moments (EDM) of the neutron 
and electron provide constraints on the phases and 9a- In the constrained minimal 
sypersymmetric standard model (CMSSM) in general the bounds on 9^ are stronger 
than those to 9a-, where 9^ ranges from O(10~^) to O(10~^) [20]. In the minimal 
SUGRA model 9^ ~ 0(10^) while 9a can be 0(1) in the smaU 9^ region [22]. 

To which extend 9b^ and 9bs arc affected by these bounds depend on the values of 
tan 13s = Vu/ Vfi, the ratio of the vacuum expectation values of the two Higgs fields in the 
MSSM. For low-to-moderate values of tan Ps in SUGRA models 9a does not shown in 
the phase of either the matrix elements {B^\n^^^^^\Bf) and {K°\n^^^^^\K°), hence 



thus 9Ba ^-nd 9bs can be safely treated as zero in this case. This means that CP violation 
asymmetries give information just for a, (5 and 7 as defined in the SM and the BSM 



(i) S{xt) ^ S{xt) + 5S{xt) 



(ii) ri2S{xt) r]2[S{xt) + 6S{xt)eJ, 
r]BS{xt) riB[S{xt) + SS{xt)AmBj 




arg(V;*,K,,), 
arg(V;*,V,,), 



(32) 
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contribution to is aligned with the tt contribution of the SM aTg{V^*^Vts) ■ In this 
context then the analysis of the UT and CP violation phases a, j3 and 7 can be carried 
out in a very similar way as in the SM, just by taking into account the contributions 
to iy*^tb) and iytdVts) without involving additional phases to those in the MSSM 
models. For large tan (3s cases it is necessary to take into account new operators whose 
contributions become relevant for this case. These contributions can be computed for 
a given theory from the effective Hamiltonian hIq^~'^^ = ^^^-^ Z]f=i Ci{iJ,)Qi, where 
Cj are the Wilson coefficients and the operators Qi are given by 

Q.^dl^.hldi^dl Q.^dlb^d^X, Q^^d!Xdib%. (33) 

The first two operators are as in the SM, while the supersymmetric contributions to 
Q3 from which chargino contributions to C^dj) are generically complex relative to the 
SM contribution and hence can generate a new phase shift in the — mixing 

amplitude. This is quite significant because C'i{fj,) oc (^^^cosjS^ . When the EDM 
constraint in these cases is imposed 6b^ becomes very small [26], [27]. 

If we like to accurately test the possible theory for physics beyond the SM, we 
have to compute all the quantities involved in the box diagrams associated to these 
mixings. Precisely the NLO QCD parameters in the supersymmetric MFV (MFSV) 
context required for Am^^, Am^^ and have been calculated under the following 
simplifications [31]: (i) the s-quark fiavour mixing matrix which diagonalizes the s- 
quark mass matrix is approximately the same as the Vckm, apart from the left-right 
mixing of the top quarks; (ii) the first and second generation s-quarks with the same 
gauge quantum numbers remain highly degenerate in masses but the third generation 
(especially t) can be significantly lighter due to RGE of the top Yukawa couplings; 
(iii) the phases 6^,3 can be safely set to zero in the entire tan f3s space once the EDM 
constraint is imposed (considered first in [25]). Then AniB^ can be expressed as 

Ams, = -^^msjlBB^VB [AsM{n) + Ah±{B) + A^^{B) + A~,{B)] , (34) 

where Asm{B) = S{xt)\V:^Vt,\' . 

The expressions for Au±{B)^ A-^±{B) and Ag{B) are obtained from the SUSY box 
diagrams. Here, H^^ xf-, ia and di represent, respectively, the charged Higgs, chargino, 
s-top and s-down-type s-quarks. The contribution of the intermediate states involving 
neutralinos is small and usually neglected. The expressions for Ah±{B), A^±{B) and 
A-g{B) are given exphcitly in -[28], [21], [15]-. 

Similarly in this context |e;^| can be expressed by: 

[e^l = C,Bk [Im Asm{K) + Im Ah±{K) + Im A^±{K) + Im A~g{K)\ , (35) 

where, the constant Cg and Im Asm are as given in Eq. (9) and the expressions for 
Im Ah±{K), Im A^±{B) and Im A~g{B) can be found in -[28], [21], [15]-. 

In the MFSV context [29] , apart from the SM degrees of freedom, only charged Higgs 
fields, charginos and a light s-top (assumed to be right-handed) contribute, with all 
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other supersymmetric particles integrated out. This scenario is effectively implemented 
in a class of SUGRA models (both minimal and non-minimal) and gauge-mediated 
models [30] in which the first two s-quark generations are heavy and the contribution 
from the intermediate gluino - s-quark states is small [22, 23, 19, 24, 21]. 

The phenomenological profiles of the unitary triangle and CP phases for the SM 
and this class of supersymmetric models can thus be meaningfully compared. Given 
the high precision on the phases a, /? and 7 expected from experiments at 5-factories 
and hadron colliders, a quantitative comparison of this kind could provide means of 
discriminating between the SM and this class of MSSM's. 

In the next subsections we follow the work of Krauss and Soff [31] in order to 
identify the NLO QCD supersymmetric contributions to the AF — 2 processes. 



7.1.2 NLO QCD-corrected effective Hamiltonian for AB = 2 

The NLO QGD-corrected effective Hamiltonian for AB = 2 transitions in the minimal 
flavour violation SUSY framework can be expressed as follows at the bosonic scale 
/x = 0(m|J[31]: 

(AB=2) _ 



eff 



47r2 



(36) 



where Qi is given in Eq. (33), 5 is a sum of the Inami-Lim functions for the different 
internal particles, S — S{xw-i xh) + S{xi, ya) ^, where 



Xw,H 



mt 



Ml 



Xi 



Ml 



(37) 



Qi(yu) and 1]2{B, fi) depend on the scale but once they are run down to the electroweak 
scale /i = Ml and the matching conditions are performed, then the NLO QGD cor- 
rection factor 772 (-B) can be expressed as [31]: 



r]2{B) = asimw) 



1 + 



4:71 



(38) 



in which n/ is the number of active quark flavours (here Uf = 5), Znj 



2A 



(0) 



2A 



(0)^ 

■I 



and 7*-°'' and (3^^ are the anomalous dimension and beta functions of QGD ^ with N(. = ?> 
and Cp = 4/3. The function D comes about from identifying the QGD correction factor 
once the Wilson coefficient is computed (in the naive dimensional regularization scheme 
using MS) and is a function of Eq. (37), D — D{xw, xh, x^^) -\- D{{xi, x^^). The 

"^Sixw, xh) = Sww{xw) + 2 Swh{xw, xh) + Shh{xh) and 

S{{xi,ya})=^ J2 E Kij^abS{xi,Xj,ya,yb) 

ij'=l a,b=l 



7 



(0) 



6 



:j(0) _ llJVc-2ra/ 



Puf — 



2G/n/ 



lUf — 



_21 -I- -i!- — -tii AT 



I"/ 



27 



factor [1 + "^^^^^-^ + Znj^j] is quite stable against variations of the functions D and 
S and has been found to be [31] and [25], around 0.89 and hence f]2{B) it is found to 
be r]2{B) = 0.51 in the MS scheme. 

The Hamiltonian given above for — mixing leads to the mass difference 

Am^, = ^{V,aV,l)%{B) S flBs, , (39) 

and analogously for Am^^. Since in the MFV supersymmetric context, the QCD cor- 
rection factors are identical for Am^^ and Attib^ then Auib^ and Ams^ have the same 
enhancement, with respect to their values in the SM. Hence the ratio AmsJ Ams^ is 
the same as in the SM. 

7.1.3 NLO QCD-corrected effective Hamiltonian for A5 = 2 

The NLO QCD-corrected Hamiltonian for AS* = 2 transitions in the MFSV framework 
has also been obtained in Ref. [31]. From this, the result for [e^l can be written as: 

|e;,| = C,BKA^X%-r]^x, + A^A" (1 - p) miK) S + rj^Six,, Xt)], 

where the NLO QCD correction factor is [31]: 



(40) 



again, as it happens with r]2{K), the factor involving the functions Zj, D and is quite 
stable and has been calculated [31, 25] to give 0.84, consequently r)2{K) is estimated 
to give rj2{K) = 0.53 in the MS scheme. 

7.1.4 Effects of MVSV for Am^^ and ek 

Here wc summarize the effects of MVSV for the transitions AB = 2 and A^* = 2 which 
with respect to the SM may be characterized by the shifts 

r]BS{xt) 772(5) S, 

V2S{xt) ^ m{K)S, (41) 

where 772 (-B), r]2{K) and S are functions of (m^±,mf2,mi7±,tan/35). In Section (7.3.2) 
we are going to use the approximations of [25] and [31] for 772 (-B) and r]2{K) 

7]2{B) = 0.51, 

il2{K) = 0.53 (42) 
and compute explicitly S for different supersymmetric inputs. 
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7.2 Scenario II 



In this context we can test models for which the BSM contribution to ex it is consid- 
erably larger than the corresponding contribution to ^niB^ and Auib^ or vice- versa. 
Hence we can expect variations of the top-W box diagram function, S{xt), different for 
i^ruB^ (or AmsJ and for ck- 

An example of this possibility occurs in models with a horizontal symmetry where 
often the CKM matrix is no longer the only source of CP violation because there are off 
diagonal elements in the soft squared matrices as well as in the trilinear terms. In this 
respect these models depart from the condition of MFV since the dynamics of flavour 
violation is no longer completely determined from the structure of the SM Yukawa 
couplings. However if the models behave in the low energy limit as a supersymmetric 
extension of the SM, i.e. all fields other than the SM ones and their supersymmetric 
partners have been integrated out at the scales where the supersymmetry is broken; 
then we can describe Am^^ and eK as in Eq. (34) and Eq. (35) respectively where we 
can have sizable contributions for example from the gauginos. 

In [45] we have constructed and analyzed the predictions of a supersymmetric model 
with an underlying S'f/(3) horizontal symmetry and with spontaneous CP violation. 
A full numerical analysis of this type of models, including the exact diagonalization of 
the soft mass matrices, needs to be performed, but since in the context of the Mass 
Insertion approximations (MI) ^ the AF = 2 Hamiltonian can be expressed [42] in 
terms of 5^^, and 5^^, we have compared the experimental bounds with the 
predictions of the model. The model puts the following hmits 



where $i is a phase of the model, preferred to be small. These hmits can be compared 
with the experimental limits, obtained in the kaon sector, of A/Im(5^^)i2 < 3.2 x 10~^ 



and Y [Ini((5^j^)i2((5^i:,)i2] < 2.2 x 10~^, which correspond to average s-quark masses 

of ihq — 500 GeV. These can be scaled as mqGeV/500, so we can see that for large 
values of rrig, we can have a large contribution to e^- However as we will see for the 
current results on possible contributions of non SM physics to e/^, these are rather 
small, limiting the possibilities of a big contribution to eK from these kind of theories. 
In the B sector the contributions can be negligible, so this is an example in which ek 
can receive big contributions, while AmB^ is left practically unchanged. 

^Defined by = ^ ^ , where A,B = L,R refers to the chirahty of the super-partners of 
the internal hnes of the relevant AF processes and are the matrices diagonalizing the Yukawa 




(43) 




couplings, e.g. Yi = V^^Y^^V^ 
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7.3 Results 



7.3.1 Procedure to extract information on bounds BSM 

The procedure to obtain the pdfs for r^^, and r^^ follows from the very simple 
observation that we can write these parameters as 

- - (44) 



where q/\mBa, qi,CK and q2,€K ^i-re functions of the Wolfenstein parameters (and of the 
correspondent AB = 2 and A^* = 2 transition amplitudes). We can extract the pdf 
for tb^ and r^^ as p{r) oc / p{p,f])\dp/dr\dfi, exactly as in the case of the parameters 
of the Classic Fit, Eq. (17). The most conservative way to obtain r^^, Ob^ and ex is 
to use only the constraint |\46|/|V^fe| and experimental constraints on a and 7, such as 
those associated with the decays B pp. pn and pn for a and B DK for 7, which 
are AF = 1 processes. However we use the constraint AtubJ Atubj^ for two reasons: 
the first one is to study the impact of the Auib^ measurement and the second one to 
probe the models for which this ratio is left invariant, such as the MFV case. In order 
to present the impact on the parameters measuring processes BSM we perform the 
following sets of fits. 

(i) Constraints iKtbl/lKbL sin2/3, AtubJ Auib^ and Atjib^ to obtain rg^, 

(ii) Constraints AuibJ Attibj and ex to obtain rAms^ ^^"^ ^b^j 

(iii) Constraints |Ki&|/|K;6|, sin2/3, AtubJ AmB^ to obtain r^^, '"Ams^ and Ob^- 

In case (iii), should be consistent with zero at some CL since we are imposing the 
sin 2/3 constraint. Prom these fits we can extract then information for bounds on the 
different examples presented in the previous section. 

7.3.2 Outputs 

(i) Results using the constraints iVubl/jVcbl, sin 2/3, AniBs/AniBj and AniBj 

The results of this fit are relevant to probe the models for which Am^^ is expected to 
be left invariant while a change in ex may happen, as it is the case of the SU (3) fiavour 
models [45], as pointed out in Section (7.2). First wc can check that the experimental 
value of ex lies within a third of the 68% CL of the output value of ex- 2.56+o!26 

xlO"-^. 

Then we see that r^^ is compatible with unity: 

r,, = 0.906lg:ii^ (45) 
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Figure 16: One dimensional probabilities for (i) ex, (ii) r^^^ and (iii) 6S{xt)eK when ex is 
not imposed constraint. 



although the relative uncertainty in this parameter is about 13%. Taking into account 
that the relative uncertainty in the hadronic parameter is of about 12%, since 
Bk = 0.707 ± 0.0066 ± 0.072, we see then that the uncertainty in the determination of 
r^^ is due to the uncertainty on Bk, since indeed we can only determine with a good 
accuracy r^^^Bx- However it is important to remark that the value of r^^^ does not 
cluster exaclty around 1, which leaves interesting possibilites of processes BSM in this 
case. For the example of SU{3) model presented in Section (7.2) we see that if this 
model is to satify the present determination of r^^^ then the s-quark masses should not 
exceed 500 GeV and the phase $i should be quite small. 

In the general context of MFV D'Ambrosio et al. [16] have classified the operators 
for which we can expect contributions BSM. In particular for the AF = 2 operators. 
On the effective Hamiltonian can be written as 



^cff = -J^ E ^nOn + h.C. ^ ^^^2 a ^tlVtj E CnQn + h.C, (46) 

J\ n 2v 27r sm Ow n 



where A is the scale at which the operators become relevant and a its effective coef- 
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ficient. The term after the arrow is the standard notation for the Hamiltonian using 
the Wilson coefficients C„ and the QCD operators Q„. The Hamihonian in Eq. (46) 
can be normahzed such that the electroweak contribution is of order one, i.e. for the 
SM asM = 1 and the scale Asm = Aq can be defined such that Ao = yt sin^ OwMw jo- ~ 
2.4 TeV. Hence we consider the ratios tb^ and r^^^, defined through Eq. (12), getting a 
measure of (Aq/A)^. In particular 

5S{x,) = ^{^)\ (47) 



Asm V A 

such that we can put bounds on the scale A at which a process BSM could take place, 
by considering a = 1. For this case, taking into account the full 68% CL range, the 
bounds for A are A > 5.84 TeV for a negative 5S{xt) contribution with respect to the 
SM contribution and A > 17.52 TeV for a positive contribution. 

(ii) Results using the constraints |Vub|/|Vcb|, AniBs/AniBj and ex 




Figure 17: The two dimensional probability (i) for r^d ^-nd 6b^ and the one dimensional 
probability (ii) for rsd- As we can see this solution agrees remarkably well with the SM, for 
which VBd = 1 and 9b^ = 0, although there could be a phase different from zero, giving new 
contributions to the phase /3 and a. 

Here it is important to compare the output value of the fit for Aitlbj to its ex- 
perimental counterpart and to the output of the complete fit to check that this fit is 
consistent with the experimental value. The output values of Am^^ for this fit and the 
complete fit are respectively 

AniB^ = (0.488 - 0.294 + 0.291)ps"\ Am^, = (0.476 - 0.196 + 0.183)ps-\ 

which are to be compared to the experimental measurement Am^^ = (0.506±0.005) ps~ 
We can see that the output value of Am^^ agrees remarkably well with its experimental 
counterpart, consequently the output value for r^^, 

= 0.983 - 0.108 + 0.091, (48) 
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is still consistent with 1. This value and its relative uncertainty (10%) however have 
decreased with respect to the value used when considering the last bound on Ams^, 
which gave r^^ = 1.095 — 0.315+0.173 with a relative uncertainty of about 23%, leaving 
some possibihties for processes BSM. We can check this also by determining the phase 
Osa: as defined in Eq. (13), 



which is compatible with the SM at 68%, although we could have an interesting pos- 
sibility of having a contribution to the angles /3 and a, through the relation Eq. (13), 
while no contribution BSM to Am^^. 

(iii) Results using the constraints |Vub|/|Vcb|? sin 2/3 and AmB^/AmBd 

The results of this fit can be used to test the hypothesis of MFV violation and 
specifically the supersymmetric case, with the extra assumptions mentioned in Section 
(7.1). In these cases it is assumed that Am^^ and |e/^| will have the same change in 
S{xt), the box diagram function of the top-W boson interaction, thus in general we 
can compare r^^ and r^^ as defined in Eq. (1). We quote here the values using the last 
Am^^ bound and the current Am^^ measurement (the pdfs of these can be found in 
figures (18 ii ) and (19 ii), which show important differences. 



where the first line corresponds to the last bound on Auib^ and the second to its current 
measurement. Using the last bound on Am^^, 5S{xt)eK and SS{xt)AmB^ agreed pretty 
well with the MFV scenario: 

SS{xt)e^ = 0.684 - 0.565 + 0.280, SS{xt)AmB^ = 0.587 - 0.950 + 0.846, (51) 

but the uncertainties were quite big. With the current AtuBs measurement, although 
there is overlap for the values of S{xt) extracted from ex and Am^^, 

5S{xt)eK = -0.373 - 0.375 + 0.343, 5S{xt)AmB^ = -0.082 - 0.720 + 0.359 (52) 

the MFV assumption appears to be weakened. 

According to Eq. (47) the limits on A from Eq. (51) are for ex'- A > 4.6 TeV. In 
the 68% CL range in this case 6S{xt)ei^ is a positive contribution. For AniB/- A > 8 
TeV for a negative contribution with respect to the SM and A > 4 TeV for a positive 
contribution. 

For the case of the supersymmetric MFV we can say more, since we have the 
change of Eq. (41), thus we can determine to which values of rj2{B), r]2{K) and S the 
parameters determined from the fits correspond to 



= (-3.35 - 3.06 + 4.03) n. 



(49) 



r^^ = 1.211 - 0.195 + 0.280, 
= 0.904- 0.121 + 0.110, 



TBa = 1-258 - 0.422 + 0.354, 
TBa = 0.933- 0.097 + 0.102, 



(50) 




S{xt) 




S, 



Vtt 



[S{xt) + SS{xt)e,] = 2.229 



1+0.343 
-0.354 



(53) 
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where we have used S{xt) = 2.2764 — 0.06466 + 0.0652, as given by the experimental 
information on Mw and rrit. Here we would like to mention the great improvement 
with respect to the values obtained using the bound Eq. (20) 

Saitib^ = 3.0006_,_Q g444, 

S,, = 3.2062;?:^8i9, (54) 

For SAniB^ there is an improvement in the relative uncertainty from 33% to 14% and 
for S^j^ from 31% to 16%. The values of 7]^ and rjB can be read from Table (2). The 
values of rj2{B) and rj2{K) are given in Section (7.1.1) in Eq. (42). 

We can determine values of 5*, with a given input of representative supersymmet- 
ric extra-constraints to the MFV hypothesis. In Table (1) we present some of these 
examples, which need to satisfy the conditions 

150 GeV < mi^ < 300 GeV, < m^^ < 600 GeV, 

< < 900 GeV, 
100 GeV < m^2, m^, < 400 GeV, (55) 
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Figure 19: In figure (i) we liave presented tlie contour levels of rs^ versus at 99%, 95% 
and 68%. In (ii) the one dimensional pdf of rs^ and in (iii) the one dimensional pdf of the 
solution for 9b^ compatible with the SM, whose 68%CL is (— 5.99, — 2.17)[°] with mean at 
-4.03° and another solution whose 68% CL is (43.01, 46.21) [°] about 44.64°. 



for rriq = rriu^ j ^ ^ . We can see that larger values of 5* correspond to small values of 
tanPs and a heavy supersymmetric spectra, except for the charginos, which can be 
somewhat lighter than the rest of the involved s-particles. 

All examples here presented. Table (1), are compatible with the results of Eq. (54), 
at 68% CL, specially those with tanPs = 8. This is quite interesting because the 
values using the last bound on Am^^, Eq. (54), had a limited compatibility, from the 
the examples presented in Table (1) the ones for tan /5s = 5 are in agreement with the 
determination of 5* from the observed values of r^^, S{xt) and 6S{xt)ejf at 68% CL, 
while the last solution for tan /3s = 8 was in agreement just at 99% CL. 

More constraints on the determination of the UT and a better determination of the 
QCD parameters will reduce the uncertainties in the determination of S from Eq. (53) 
and hence will help to possibly rule out some possibilities within the supersymmetric 
MFV scenario. 
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rriH 


nig 




mi 






tan ps 


b 


115 


880 


400 


150 


170 


189 


5 


2.56 ± 0.08 


lib 


880 


400 


150 


170 


189 


5 


o ro-0.08 
Z.OZ_|_o.o6 


1 lO 


ooU 


Ann 


1 "^n 

lOU 


"^nn 


OOU 


o 


9 Afi 4- n n7 
z.^o zn u.u / 


115 


700 


350 


160 


170 


189 


5 


2.69 ±0.09 


200 


700 


350 


160 


170 


189 


5 


2.63 ±0.09 


115 


700 


350 


160 


170 


189 


8 


2.58 ±0.09 


115 


700 


350 


160 


300 


380 


8 


2.35 ± 0.07 



Table 1: Examples of input parameters to determine S, all the masses are expressed in GeV. 

8 Conclusions 

The goal of this note was to study the impact of the present limit on Am^^ on the 
determination of the unitary triangle (UT) parameters within the context of the SM 
and also to put bounds for certain classes of processes BSM for which the current UT 
analysis it is useful. 

The first version of this note used the bound on Atjib^ > 14.5ps^^ and the bound of 
AniBe > 16.6ps^^ in order to begin to study the impact of Am^^ on the SM and other 
models. These bounds were set by combining all the Atub^ information available at 
the end of 2005, including the one reported by the CDF and D0 collaborations prior 
and after Tevatron run II, respectively. The first measurement of Am^^ = 17.33 ±0.4, 
which has an uncertainty of 2% its value, of course makes a stronger impact on the 
determination of the CKM parameters and possible BSM scenarios, being limited only 
by uncertainties from hadronic computations. 

Using as an input the bound Am^^ > 16.6ps~^ we obtained an output value of 
AniB^ — 19.96;^Ii;qi which agreed remarkably well with its output value with out using 
AniB^ as a constraint: Am^^ = 21.08l,l2!y3. However we can see that the measurement 
of AiTiB, = 17.33 ± 0.4 differs from this estimation, although it is not incompatible. 

We tested the consistency between the CP conserving processes - |Ki6|/| K;6|, Atub^ 
and ArriBa/ Attib^ ~ violating ones - €k and sin 2/3. With the last bound on 

AniB^ and the used set of experimental information there was an slight disagreement 
between the resulting value for sin 2/5 considering just CP conserving processes, which 
in this case renders a value of sin 

2^ CP-conscrv. ^ (0.751,0.841), and the experimental 
reported value of sin2/?^^P- = (0.623, 0.751), both intervals at 95% CL. With the mea- 
surement of Attib^ = 17.33 ± 0.4 the compatibility it is just at 99% CL since for this 
case we have sin 2/3 ^P-conserv. ^ (0.742,0.856). 

An interesting possibility is that there is a beyond the SM (BSM) contribution to 
sin 2/3 which it may be expressed as sin 2/3 = sm{2(3^^ ± ^bJ- We have checked that 
this case it is indeed possible by taking into account the results for this phase. From 
the parameterization rs^e^^^d = {Bl,\nJs''^\Bl^)/Bd,s\n^s\Bd,s) and the results of 
Eq. (49), we obtain Ob^ = (-3.35 - 3.06 ± 4.03j[°] and = 0.983 - 0.108 ± 0.091, 
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both at the 68%. CL. 

Theories satisfying the minimal flavour violating (MFV) conditions, for which 
iKtfel/lKbl s-nd ArriB^/ AniB^, used for the unitary triangle fits in the SM remain in 
the same form, can currently be tested at the same level of precision at which the 
analyses in the SM are carried out. In particular a supersymmetric version of the 
MFV scenario where the first two generations of s-quarks are taken to be heavy and 
degenerate, allowing just the contribution of the third generation, iji and II can be 
tested in more detail since the change at the QCD NLO is known [31]. We have tested 
this possibility for some samples of supersymmetric spectra, finding an agreement with 
most part of the supersymmetric examples considered, specially for tan (3s near to 5 (in 
this scenario tan cannot be large). 

The first measurement of Atub^ = 17.33 ± 0.4 and its relatively small uncertainty 
(2%) has provided a good estimation of the region of compatibility of the MFV su- 
persymmetric scenario, a better determination Ame, and |T46| and \Vcb\ will help to 
narrow the uncertainties in all the parameters related to it in the UT analysis, specifi- 
cally the parameters related to observations BSM and hence helping us to disentangle 
more possibilities in this respect. After the first version of this work was presented 
some related works, specially considering the impact of Am^^ on the MFV scenario or 
other BSM scenarios appeared in the literature [46] and the unitary triangle updates 
after the first measurement of Am^^ [47]. 

Acknowledgments 

L. Velasco-Sevilla would like to thank F. Krauss for the clarification of the notation 
of his work in hep-ph/9807238, M. Voloshin for discussions and also take here the 
opportunity to deeply thank A. Stocchi for help understanding the bayesian approach 
when developing the code for the determination of the CKM parameters. L. Velasco- 
Sevilla is supported in part by the DOE grant DE-FG02-94ER. 



37 



A Experimental Inputs 



Fixed Parameters 



P^rQiTlPtpr 

A. CHi± dlllC ud 


\/a 1 n p 


R pf prPTi pp 




i.iuuoy A iu vje V 






(80 425 ± 038)GeV 


n 


fK 


(0.1598 ±0.0015)GcV 


u 


rriK 


(0.49765 ±0.00002)GcV 


u 


ArriK 


(3.4606 ± 0.006) x lO-^^GeV 


u 


\eK\ 


(2.280 ±0.017) X 10-3 


u 


ritt 


(0.574 ± 0.004) 


[37] 


TUB, 


(5.2794 ± 0.0005)GeV 


[33] 




0.55 ±0.007 


a 




(5.375 ± 0.0024)GeV 


a 



Table 2: Input values of fixed parameters. 



Fitted Parameters 



Parameter 


Value ± Gaussian errors 


Flat errors 


Refer en. 


iKbKincl.) 


(41.6 ±0.7) X 10-3 




[32] 


V^b (excl.) 


(41.3 ± 1.0) X 10-3 


1.8 X 10-3 


a 


iKbKincl.) 


(43.9 ±2.0) X 10-^ 


2.7 X 10-^ 




iKbKexcl.) 


(38.0 ±2.7) X 10-^ 


4.7 X 10-^ 


11 


IK.I 


0.2258 ±0.0014 




a 


Bk 


0.79 ±0.04 


±0.09 


[33] 


rUc 


(1.3±0.1)GeV 






rUt 


(161.5 ±3.0)GeV 




[35]* 


Vcc 


1.38 ±0.53 




[36] 


Vet 


0.47 ±0.04 




[37] 


sin 2(3 


0.687 ±0.032 




[33] 




(0.506 ±0.005)ps-i 




[34] 


fBs\/BBs 


(0.276 ± 0.038)GeV 




[41] 




1.21 ±0.04 


±0.06 


[38] 




(17.33 ±0.4) ps-i 







Table 3: Input values of the fitted parameters. * Taken from [35] and calculated at the pole. 
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B Tables of outputs 



Outputs 



L CtiCtiliCtCi 


/\td n ^ 1 f\ f\ Tic; 


Nn Am o 


A™ _ fiy oo _)_ Q „o -1 


P 


n 1 — O.U^2 
U.ZiZ_,_o.o34 


U-^<JO+0.034 


n 900^0. u:^i 
U.ZUU+0.022 


V 


U. 000+0.028 


U-<J40_,_o.025 


U.JDO+0 015 


L^mB^ [ps J 


J-y- -1-0+1.35 


Zi.UO_,_2.73 


1 7 49^0-38 
i i .4/:_,_o.27 


\V AIWA 


u.uu»/ +0.005 


u.uu»o_|_o 005 


u.uu» / _|_o.o05 


Am^^ [ps-^] 


U.OU4_,_o 057 


o.483;|]:li 


0-476+0.183 


sin 2/3 


Q 7P-P--0.019 
U. i OO+0.023 


n yc-c- -0.024 
U. <00+o.023 


U- '00+0.028 




98.90;t^9 


i02.98;?:g 


97.49;|i[''] 


in 


59.39;i:85 


55.32_|_5;3g 


60.97;t°^ 


Rb 


n 41 Q-0.008 
U.4iy_^0.013 


'-'■^-^-'-+0.020 


o.4i6;Hi^ 


Rt 


U.OOO_,_o.o36 


0.836+o:o35 


n oiy^— 0.022 
O-oi 1 +0.022 


Bk 




0.786;o:o47 


n 7c;c;-0.041 
U. / OO+0.042 


fB.^BsAGeV] 


0.226+o!oi4 


U.ZO4+0.016 


n 99/1-0.016 
U.ZZ4+0.014 



Table 4: Output values of the fitted parameters. The three columns correspond respec- 
tively to the complete Classic Fit using the bound Ams^ > 16.6 ps without using any 
information on Auib, and finally using the measurement of Auib, = (17.33 ± 0.4) ps 



Outputs 



Parameter 


(i) 


(ii) 


(iii) 


P 


U.iOO_,_o.030 


o.209;H^° 


0.184J:°^^ 


V 


n Q7C,-0.020 
U.O/0_|_o.oi8 


o.379;H{I 


A o7p:-0.018 

u.o ( 0+0.017 


Amsjps-^] 


i7.37;°:i 


17.48;°i? 


1 7 "^fi-o-^s 

1 / ■O0_|_o.38 


\Vuh/V,b\ 


n 0070002 

O-Us" 0.003 


o.ioi^gi 


0.0970.0020.003 


Ams^ps"^] 


U.OUO_|_o.i58 


o.488;°:if 


0.506;°:1I 




(2.56;[J:i) X 10-3 


(2.30;H^) X 10-^ 


(2.57;°:i) X 10-3 


sin 2/3 


0.766;H1I 


o.789;Hl 


o.766;°:[ii| 


«[1 


94.79;t^? 

63.6i;li? 


97.36;|«? 


94.98;|:in 


in 


60.99;|27 


63.65;|:|;^ 


Rb 


U.4ZU_,_o.o21 


0.434;H1^ 


o.42o;°:°i^ 


Rt 


n Q7Q— 0.054 
O-O' <J+0.048 


o.876;H?^ 


0.897;^! 


Bk 


0.707;°:°?^ 


A 740O.O36 

"J- ' ^"^0.030 


0.706;°:°?° 


fB.^BBAGeV] 


0.219;°:°!^ 


A 99t--0.008 
U.ZZd_|_o.oo7 


0.219;°:°}^ 



Table 5: Output values of the fitted parameters. The three columns correspond to the three 
cases presented in Section (7.3.2): (i)Using the constraints |V^6|/|Fcfe|, sin2/?, AtubJ Atub^ 
and AmB^; (ii) iKbl/lKbl, AtubJ Ams^ and ek and (iii) iKbl/lKbl, sin 2/3 and AmBjAmB^- 
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